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We revisit the topic of triple-product asymmetries which probe CP violation through differential 
distributions. We construct distributions with well-defined discrete symmetry properties and char¬ 
acterize the asymmetries formed upon them. It is stressed that the simplest asymmetries may not 
be optimal. We explore systematic generalizations having limited reliance on the process dynam¬ 
ics and phase-space parametrization. They exploit larger fractions of the information contained 
in differential distributions and may lead to increased sensitivities to CP violation. Our detailed 
treatment of the case of spinless four-body decays paves the way for further experimental studies. 


I. INTRODUCTION 

The fully differential rates of some multibody meson 
decays are being more and more accurately measured. 
In the search for new sources of GP violation, such pro¬ 
cesses present several advantages. They often feature a 
rich variety of interfering contributions from which dif¬ 
ferences in GP-violating— weak —phases could manifest 
themselves. In addition, the multiplication of measur¬ 
able independent four-vectors permits the construction 
of so-called triple-product observables. These have a cou¬ 
ple of interesting characteristics. Unlike total rate asym¬ 
metries between GP-conjugate processes, their sensitivity 
to small differences in GP-violating phases is not condi¬ 
tioned by the presence of GP-conserving— strong or uni¬ 
tary —phase differences. They can also be measured us¬ 
ing untagged samples in which GP-conjugate processes 
need not be distinguished, provided their fractions are 
equal. 

In this paper, we explore the variety of possible triple¬ 
product observables. The ever-increasing amount of data 
collected allows finer details of the differential distribu¬ 
tions for which they are proxies to become measurable. 
We stress that the most common asymmetries may not 
be the most sensitive ones, due to cancellations in phase- 
space integrals. As much as possible, we would like to 
abstract our treatment from the particular dynamics of 
the studied process. In many multibody decays, only 
phenomenological descriptions of various degrees of ac¬ 
curacy are achieved. They may not capture all the fine 
details of interfering contributions which could reveal GP 
violation. A systematic procedure that is less likely to 
miss unpredicted forms of GP violation is therefore de¬ 
sirable. Although we will mostly focus, for concreteness, 
on four-body meson decays involving spinless particles, 
our discussion has a wider range of application. 


A. Differential CP violation 

Let us consider two transitions of amplitudes 
M{{Xi,Pi}) and M{{Xi,Pi}). They involve an equal 


number particles respectively labeled by i and i, with he- 
licities Xi^^ and four-momenta We would like to per¬ 
form a comparison of these two amplitudes phase-space 
point by phase-space point so we take Aj = Xi as well as 
Pi = Pi- 

If these two processes are GP conjugate of each other, 
with i = GP[i], GP violation at any phase-space point 
takes the form of a difference between the squared moduli 
of 

Af({Ai,pJ) and M{{Xi,Pi}) 

where p = P[p] is the parity conjugate of the momentum 
p. Testing GP conservation phase-space point by phase- 
space point thus implies a comparison of the differential 
rates of two processes involving GP-conjugate particles 
of identical helicities but opposite three-momenta. 

It reveals useful to define an operator, called motion re¬ 
versal and denoted here by T, that reverts both momen¬ 
tum and spin three-vectors [1, 2]. Its action on helicities 
and momenta is thus identical to that of GP and it can 
be viewed as the unitary component of the antiunitary 
time-reversal operator T. It is therefore sometimes called 
naive T. In general, the amplitudes above can then be 
decomposed into two pieces that are respectively T-even 
and T-odd [3]: 

M{{X^,Pi}) = Me{{X„Pi}) -k Mo{{Xi,Pi}), 

M{{X^,pi}) = Me{{X^,pi}) -k Mo{{Xi,pi}) 

— AIe({A2, Pz}) Ado({Az, Pz}). 

Those two terms can receive several contributions whose 
absorptive parts [4, 5] take the form of GP-even phases 
S. One can then write 

Ade({Az,Pz}) = ai 

A^o({Az,Pz}) =aS C(^o+[‘P4+r/2])^ 

Ado({Az,Pz}) = a^ 

with implicit summation over the j, k indices, and real 
®e’o> ^e’oi Pe’o functions of the helicities and momenta 
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The above conventions imply that all CP viola¬ 
tion is encoded in the CP-odd phases (pi’ll. When they 
vanish, 

({-^2 ) Pi}) — T Ade ({ , p^}) , 

•Ado ({ Aj, Pi }) — -Ado({ Ai, Pi }) , 

SO that the CP-conjugate rates are identical, phase-space 
point by phase-space point. As the physical amplitude is 
defined up to an overall phase, a departure from zero for 
differences in these is what we are after. 


nonvanishing differences in CP-even phases S while the 
T-odd-CP-odd does not. This can be understood as, in 
the absence of absorptive part to the amplitude, T is 
equivalent to T so that CPT conservation imposes any 
CP-odd quantity to also be T odd [7]. 

On the other hand, the T-odd-CP-even piece of the 
differential rate could be used to isolate relatively small 
differences in CP-even phases S, in the absence of CP- 
odd phase p. It can thus serve to better understand 
final-state interactions. 


B. CP violation without CP-even phases 


C. Untagged samples 


The T-transformed differential rates are obviously ac¬ 
cessible experimentally since the measured momenta can 
be artificially reversed. For processes involving only 
scalars in their initial and final states, T is actually equiv¬ 
alent to parity conjugation P. The measured differential 
rates of any pair of CP-conjugate processes can there¬ 
fore be decomposed into four pieces of definite T and CP 
transformation properties: 


Another remarkable characteristic of the T-odd-CP- 
odd part of the differential rate is that it can be measured 
with samples which contain an equal number of events 
from CP-conjugated processes. It can also be evaluated 
in the decay of self-conjugate states like the Z and h 
bosons, or any Majorana fermion. This can be under¬ 
stood by rewriting the T-odd-CP-odd differential rate 
defined in Eq. (2) as 


dF 

even 

-•-■odd 

I±T 

I±CP 

dF 

(2) 

I-t /I-kCPT 

dFA 

dd) 

even 

'~"^"odd 

2 

2 

dd> 

2 2 

d$y 


with the shorthand = {Xi,pi}. 

For simplicity, let us assume there are respectively two 
and one contribution(s) to the T-even and T-odd parts 
of the amplitude in the process under scrutiny: 


M{{X^,Pi\) = al + ial 

M{{\,pi}) = al 




All functions of the phase space are evaluated at {Xi,pi\. 
Note the convention of Eq. (1) causes the appearance of 
a factor of i in front of the T-odd term. Up to a flux 
factor, the squared modulus of this expression and of 
its CP conjugate provides us with the differential rates 
which can be decomposed as prescribed in Eq. (2): 


dd) 

dd) 

dd) 


T-even 

CP-even 

t-odd 

CP-even 

t-even 

CP-odd 


oc alai + al al + al al 

+2al al cos{5l - 51) cos{pI 

oc 2alal sin((5} - cos((^} 

+2al al sin((52 -(ji) cos{pl 

oc -2 a} al sin((5} - 5l) sin((p} 




using the fact that T is an involution: T = 1. It only 
involves d(F -I- Fj/dd) evaluated at the phase-space point 
{Ai,Pi} and at its T conjugate {Xi,pi}. 

Other discrete symmetry operators can be introduced. 
In particular, let us denote a permutation of the external 
particles as E{ii, i 2 ,..., = {E[ii],E[z 2 ],..., E[z„]}. 

For transitions involving a self-conjugate subset of ex¬ 
ternal particles, there is an especially relevant permu¬ 
tation E* that takes each particle in the subset to its 
CP conjugate. For example, E*{Ar+, 77“, 7r+,7r“} = 

{77“, K'^, 7r“, 7r“'"}. 

A part of the differential rate that is odd under a per¬ 
mutation E can also be used to test CP conservation with 
samples containing an equal number of events from CP- 
conjugate processes: 


I-E /I-kCPE dFA 

2 2 d$j' 

However, resorting to such samples is only desirable when 
a subset of the particles involved is self-conjugate. Ex¬ 
perimentally, the tagging that discriminates between the 
CP-conjugate processes then comes with an efficiency 
cost. Importantly, without tagging, what is then actually 
measured is 


T-odd 

oc 2 a} a} cos((5} - (i}) sin(vj} - pi) 

CP-odd „ , „ , „ , 

+2al al cos{Sl - 5}) sin((p} - pi). 

The last two expressions above vanish in the CP limit. 
There are thus two distinct kinds of CP-violating differ¬ 
ential rates [6]: the presence of the T-even one requires 


dd) 


I-kCPTE* dF 
2 d$' 

In an untagged sample, one can therefore measure two 
CP-odd differential rates that are either T-odd-E*-even 
or T-even-E*-odd: 
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I ± T I T E* I - CP dr 
^ 2 2 d$ 

_ i±t Ite* /i + cpte* dr\ 

“ ^ ^ 2 d$ J ■ 

Some asymmetries of either kind were measured by 
the LHCb Collaboration in its study of the —> 
K^K~'k^tt~ decay with an untagged sample [8] (see dis¬ 
cussion in Section II.H). 

On the contrary, the differential rates of identical T 
and E* parities are CP even in an untagged sample: 

I±t I±E* /I-kCPTE* dr\ 

^ 2 2 d $) 

_ I ± t I ± E* I -k CP dP 
“ ^ 2 2 d$’ 

As in the tagged sample case, they provide a handle on 
the CP-even phases. 


D. Integrated observables 


No phase-space integration or spin averaging is in prin¬ 
ciple required to test for the existence of CP-violating 
phases. Such procedures are only applied because of 
practical constraints like finite statistics. The total rate 
asymmetry is constructed upon the T-even-CP-odd dif¬ 
ferential rate 



d<i) 


T-even 

CP-odd 


(3) 


A second family of observables can be obtained from in¬ 
tegrals of its T-odd-CP-odd homologue 


fdd> f(d>) 


dd) 


t-odd 

CP-odd 


(4) 


with some T-odd function f (<&) without which the phase- 
space integral would vanish. Similarly, any T-even func¬ 
tion g(d)) could be inserted into the T-even-CP-odd in¬ 
tegral to construct observables sharing the properties of 
the total rate asymmetry. 

As a product of a T-odd kinematic function with a T- 
odd-CP-even differential rate, the observables of Eq. (4) 
are T even and CP odd but have not definite T transfor¬ 
mation properties. 


Dot products and antisymmetric contractions of four- 
momenta and Pauli-Lubahski spin vectors (respectively 
denoted by p and w) have definite P parities. The P-even 
combinations are: 

Pl-P2, WI-W2, 

PiPIPzwI, and 

while the P-odd ones are: 

Pi ■ W2, 

U V (7 P U V G P U V G p 

PiP 2P3P4^ ^pvpaPiP2W3Wl, W^W^W^W^. 

The sensitivities to discrete symmetry violation of ob¬ 
servables having definite P and T transformation prop¬ 
erties, in the presence or absence of absorptive parts in 
the amplitude, are listed on p. 519 of Ref. [6]. 

The completely antisymmetric Lorentz structure 
can originate directly from Lagrangian couplings like 
i^pupaF^'^, or arise in the presence of chiral fermions, 
since 7® = jr ■ Because it is completely 

antisymmetric, however, a necessary condition for the 
presence of a T-odd part in an amplitude is the 
availability of four independent and distinguishable four- 
vectors. In a process involving scalars or particles of un¬ 
measured spins, at least five external momenta are there¬ 
fore required. 

In a reference frame where = (a°,0), the com¬ 
pletely antisymmetric combination of four four-vectors 
C/ii/ptr b'' d 'reduces to a a° 6 • (c x d) scalar triple 
product (for eoi 23 = +!)• The observables constructed 
from the T-odd parts of the differential rate are therefore 
customarily called triple-product asymmetries. A signifi¬ 
cant amount of effort, both theoretical and experimental 
has been devoted to their study. A triple-product asym¬ 
metry has been measured in AT/, —>■ 7r'''7r“e“'“e“ [9] and 
applications are also found in heavy-meson [10-26] [27- 
31], baryon [32, 33], top [34], Z [35], Higgs [36-39], and 
beyond-the-standard-model [40] physics. 


F. Asymmetries 


The simplest up-down triple-product asymmetries are 
based on the sign of one of the constructible triple product 
(see Eq. (4)) 


f(4>) = sign{ Cp^pcr a^b’^dd'^ }. 
The usual quantities defined in the literature 


(5) 


E. T oddity and triple products 

There are two tensors available to construct Lorentz in¬ 
variants from spin and momenta four-vectors. The met¬ 
ric Ppu leads to T-even contractions like invariant masses, 
and the completely antisymmetric e^^po- produces T-odd 
combinations of four-vectors. 



/d4> ^(4>) [ 

rd$ [ dr I T-even 
1 d$ IcP-eve 


dr I T-odd 

d$ I CP-odd J 

I T-even -i 

I CP-odd J 


are ratios of integrated T-odd and T-even differential 
rates and have no definite CP transformation properties. 
The converse could only be argued when differences of 
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CP-even phases are proven vanishing. In the notations 
of Section LB, 


« 2 ai a’: sin [{Si - - ^D] 


then actually becomes a probe for small differences in 
CP-odd phases On the contrary, 


is always CP odd. Ap is occasionally defined as the 
CPT conjugate of Aj. and has then a sign opposite to 
Af = CP[^rf,] defined here. With this alternative con¬ 
vention, becomes a sum. Other asymmetries were 
for instance listed in Ref. [26]. Instead of ^1?^, one may 
consider 



/d4> m 


dr I T'Otici 
d¥ I CP-odd 


/dch 


I T-even 
d^» I CP-even 


This choice corresponds to the more common one when 
the total rate asymmetry of Eq. (3) vanishes. Using 

the T-even-CP-odd and T-odd-CP-odd families of ob¬ 
servables can be kept independent. Uncertainties in the 
relative abundance of the two CP-conjugate initial states 
can however make the use of -4?^ experimentally prefer¬ 
able. 


To avoid dilutions in the integral of Eq. (4), the functions 
chosen should ideally change sign wherever the T-odd- 
CP-odd piece of the differential decay rate itself changes 
sign. The bins’ boundaries should also be placed there. 

The question of what set of f (4)) functions would yield 
the best sensitivity to CP violation is nontrivial and de¬ 
pends on the process at hand. Actually, when the form of 
the differential decay rate is known with confidence, one 
may rely on an unbinned likelihood fit to the data for ex¬ 
tracting CP-violating parameters. Such amplitude analy¬ 
ses have notably been carried out for several R-meson de¬ 
cays: e.g., for —>• K~K~, dominated by a (/)0 

intermediate state [29], or for —>■ K~^K~, dom¬ 
inated by a (j)K*^ resonant intermediate state [27, 31]. 

Trustworthy parametrizations also make it possible to 
determine the asymmetries relevant in the study of the 
CP-odd phases that might appear in perturbative pro¬ 
cesses like h or e+e“ —)■ [39, 44- 

48]. Observables of optimal statistical significance can 
then also be determined [49]. 

In the hadronic decays of heavy mesons however, the 
parametrization provided by a resonance model is only 
phenomenological and, although it may capture accu¬ 
rately enough the main features of the studied process, 
new sources of CP violation may only be observable 
in finer details. Using tests of CP violation that have 
a limited reliance upon the process dynamics and its 
parametrization is therefore desirable. 


G. Dilutions and f(4>) sets 

The ‘sign’ function used in Eq. (5) is not the only pos¬ 
sible weight function 7(4)) that could be used in phase- 
space integrals of T-odd differential rates. This choice 
amounts, experimentally, to counting events in regions of 
phase space. Moreover, the adjunction of any T-even fac¬ 
tor in the ‘sign’ argument besides the antisymmetric con¬ 
traction would obviously yield other po¬ 

tentially interesting observables. Using a basis of T-odd 
functions on 4), it is also possible to decompose the T- 
odd-CP-odd differential rate in moments (see Refs. [41- 
43] about the method of moments). As in Ref. [28], 
a binning of the phase space could also be defined and 
a chi-squared test carried out to assess local departures 
from zero in the T-odd-CP-odd piece of the differential 
rate. This would correspond to choosing, for the 7(4))s, 
a set of characteristic functions that evaluate to 1 in one 
bin and vanish elsewhere. At least three categories of 
7(4)) functions can thus be used to describe the T-odd- 
CP-odd piece of the differential decay rate: 

— ‘sign’ functions defining a signed partition of the phase 
space, 

— a T-odd basis on 4* providing a decomposition in mo¬ 
ments, 

— characteristic functions defining a phase-space bin¬ 
ning. 


II. SPINLESS FOUR-BODY DECAYS 

Four-body decays involving only spinless particles are 
simple examples of processes in which four independent 
four-vectors can be measured. In these cases, T is equiv¬ 
alent to P and there is actually one single independent 
antisymmetric e^j/po- contraction which involves the ex¬ 
ternal particles’ four-momenta. All T-odd functions of 
the phase space are built upon it. 

In the following, we will focus on this simple case and 
investigate how to define appropriate signed partitions 
(or binnings) of the phase space. For concreteness, we 
will often refer to the specific de¬ 

cay. Its differential rate, as well as the one of the corre¬ 
sponding CP-conjugate process, has recently been mea¬ 
sured with an impressive accuracy by the LHCb Collab¬ 
oration [28]. Note we will only consider time-integrated 
quantities while the LHCb Collaboration also recorded 
the time dependence of the decay rate. 

A. Phase-space parametrization 

Hadronic multibody decays often receive contributions 
of various topologies. The ones so far measured in 
—>■ K~Tr^Tr~ are displayed in Table 1. A given res¬ 

onance structure would be most appropriately described 
with a parametrization of the phase space that includes 
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Intermediate states 

Br xl0'‘ 

ipp°)s, p^K^K~, ->• 7r+7r“ 

{Pp°)d, P^K^K~, p° —^ tt+tt” 

{K*°K*°)s, K*° 

9.3 ± 1.2 

0.83 ±0.23 

1.48 ±0.30 

P{tv'^tt~)s, P —>■ K'^ K~ 

2.50 ±0.33 

(A'“7r+)p(A'+7r“)s 

2.6 ±0.5 

K+K-, K+^ K*°'k+ 

K-K+, K- K*°tv- 

1.8 ±0.5 

0.22 ±0.12 

K+K-, K+^p°K+ 

KpK+, Kjj 

1.14 ±0.26 

1.46 ±0.25 

A:*(1410)+A-, A'*(1410)+ K*°tt+ 
A:*(1410)"A:+, A^UIO)" 

1.02 ±0.26 

1.14 ±0.25 


TABLE I. The different measured contributions [50] to 
Br(D° K+K-tv+tv-) = (24.3 ± 1.2) x lO"'^ as listed 
by the Particle Data Group [51], the corresponding decay 
topologies—or resonance structures—, and branching frac¬ 
tions. The S, P, and D indices indicate the partial waves 
in which the particle pairs are produced. 



The standard Cabibbo-Maksymowicz parametriza- 
tion [52] of the phase space can be adopted to describe 
a four-body decay of 0 —)■ ab ^ (12) (3 4) topology. It 
is based on two invariant masses and (which be¬ 
come constants in the narrow-width approximation) and 
three angles (see Fig. 1). In the a and b subsystems’ rest 
frames, the orientations of the final-state particles’ mo¬ 
menta are respectively characterized by 9a and 9b, com¬ 
prised in the [0, tt] interval. The relative orientation of 
the planes formed by the two pairs of momenta is mea¬ 
sured by ^ G [—tTjTt]. Note that 9a and 9b are T-even 
while (j) is T-odd (and P-odd). The whole (j) dependence 
of a differential distribution of definite T transformation 
properties can thus be obtained from the [0, tt] interval. 
The angle (j) also determines the sign of the triple product: 

P1P2P3P4 = ^mamb^J\{ml, ml, ml) s 9 a s 9 b scj). 

where A(x^, y^, z^) = {x + y + z){x + y — z){x — y + z){x — 
y — z) is the usual Kallen function, and mi_ 2 , 3,4 have 
been neglected. We will occasionally use shorthands like 
C(j) =cos(j), s^9 = sin^ 0, s29 = sin(20). 


B. Differential decay rates 

For a decay to four spinless particles forming two in¬ 
termediates states of angular momentum ja, jb, the am¬ 
plitude can be expressed in terms of spherical harmonics. 
With a spinless initial state, the two intermediate states 
have equal helicities A. We can therefore write 


FIG. 1. Parametrization of the phase space ofaO—1234 
four-body decay privileging the (12) and (3 4) subsystems. 
The momenta of the final-state particle pairs are pictured in 
their joint rest frames. 


the invariant masses in which resonances occur. Such a 
description would likely be the most sensitive to the inter¬ 
ferences between the several partial-wave contributions 
to that topology. On the contrary, the effects of the res¬ 
onances occurring in other invariant masses would be di¬ 
luted. Therefore, at this point already, the parametriza¬ 
tion of the four-body phase space challenges our aim at 
a description independent of the process dynamics. 

One may first consider the partial-wave decomposition 
of the dominant resonant intermediate state. Let us here 
focus on the 0 —?> o 6 —> (1 2) (3 4) topology found in the 
(pp^ resonant contribution that accounts for about 40% 
of the —> K~^K~TT~^Tr~ branching fraction. Repeating 

the analysis that follows for different parametrizations 
would be required to obtain better sensitivities to other 
decay topologies. A (14) (2 3) pairing would for instance 
allow to better probe CP violation involving a K*^K 
resonant intermediate state in the 
decay. 


M=47r Ai-^^iml,ml)Y>li9a,P)Y^^^i9b,0r 


with |A| < mm{ja, jb), and partial-wave amplitudes 
mass dimension —1. General expressions for 
n-body phase spaces and arbitrary spins can be derived 
from Refs. [53-55]. Our normalization is chosen such 
that the squared amplitude integrated over the 0a,fc and 
(p angles takes the form 


r dc9a dc9b d(p |^|2 
J 2 2 2tt 




In the ja = 1 = jb case relevant for the pp^ intermedi¬ 
ate state of —?> K'^K~ 'TT+TT , one can define the linear 
polarization amplitudes 

Ao^Al’\ 

where, for conciseness, we omitted the j, dependences. 
The amplitude can then be written as 


1 ^11 A± 

-M = Aq c9a c9b + s9a s9b cp s9a s9b sp 
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where the last term is T odd because of its sip depen¬ 
dence. Its factor of i respects the phase conventions of 
Eq. (1). Denoting by ^o,||,-L the linear polarization am¬ 
plitudes of the CP-conjugate process, the corresponding 
differential decay rate can be decomposed, as described 
before, into four pieces of definite T and CP transforma¬ 
tion properties: 


2 mo dP 
^ d$ 


T-even 

CP-even 


2 mo 


2mo 


dP 

d<i) 


dP 

dd) 


t-odd 

CP-even 

T-even 

CP-odd 


2mo 

“ 9 “ 


dd) 


t-odd 

CP-odd 






RejtloAJi -|- Aotl|J} 

ivi 

lm{A±AQ + A±Aq} 

iv^ 

Im{Aj_An -|- A_L^n} 


s20a s26t c4> , 
s29a s20{, scj) 
S^9a S^9b s2(j), 






C 9 a C 9 b 


S^9a S^9b c2(j) 


|Axp-|ix| o. 2 


S 9a S 9b s2(f) 


Rel^oAn — yfottn} 


4 V 2 

ImjAxAg - ixig} 
4^2 


- s29a s29b c(j ), 
s20a s29b S(j) 


lm{A±Al - A^An 
-L s^Ob S2(j), 


where 


^ VA(m^,mg,m^) A/A(mg,mj,m^) 

87rmo 87rm^ 

^A(m^, m|, ml) dc0a dc0h dcj) dm^ dm^ 
STrrn^ 2 2 27r 27r 27r 

The linear polarization amplitudes may receive different 
contributions each having a CP-even phase and a CP- 
odd phase ifi\: 


A (^2 ^2\ _ \ ^ i (^2 ^2\ i\5l^(m‘^ 

Ax(m„,mj = 2^ax(m„,mj,) e ‘ 


= ± ^ -(<53f - 4) “n (<^x - 

lm{AxA^±AxAp}/2 

= +^aVv co"(<5^-^^v) srn(<Px-‘^i^)- 

hi 

In this specific example, again, the different pieces of the 
partial rate exhibit the sensitivities to the CP-even and 
-odd phases described earlier. 


C. Beyond the most common observables 


Interestingly, with a single resonant intermediate state 
having ja = 1 = jb, the total rate asymmetry based on 
the integral of Eq. (3) vanishes when the Aq coefficient 
receives contributions of identical phases, or one single 
contribution. (The terms involving other linear polariza¬ 
tion amplitudes vanish upon phase-space integration.) In 
such a case, only could a differential rate study provide 
information about CP violation. 

Without assumption about the presence of identical 
phases, the most common up-down integrated asymme¬ 
try based on the sign of the triple product 


fdd) sign{s(()} 


dd) 


f-odd 

CP-odd 


also vanishes in this simple case. This illustrates—in an 
extreme way—that phase-space integration may result in 
losses of sensitivity to CP-violating phases. A nontrivial 
phase-space-dependent T-even factor in the T-odd-CP- 
odd differential rate can make it change sign where the 
triple product does not. 

Such dilutions can obviously be overcome when a trust¬ 
worthy parametrization of the differential rate is known. 
Taking seriously the simplified parametrization of the 
decay presented above, the bare examination of the 
differential rates indicates that more information about 
CP-odd and -even phases is contained in the piecewise 
integrals of Table II upon which asymmetries could be 
constructed. 

However, as already stressed, the parametrization of 
heavy mesons’ hadronic decays is only phenomenologi¬ 
cal and may miss some fine interference details that have 
the potential of revealing new sources of CP violation. 
We would therefore wish to adopt a more systematic ap¬ 
proach that does not rely on strong theoretical assump¬ 
tions about the process dynamics. 


D. A first look at the data 


for real-valued a^, ^x, ‘fx ^-nd X = 0, ||, or T. The 
corresponding CP-conjugate quantities are then given by 
Ax = J2i so that 

Re{AxA;.± AxA^}/2 


This point can be made more concrete using the re¬ 
cent experimental study of the —>■ K~^ de¬ 

cay. The LHCb Collaboration displayed in Ref. [28] the 
measured ttt-k+k-j cos9t^, cos9x and (j) distri¬ 

butions for both and as well as sin (p > 0 and 
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FIG. 2. Decomposition of the measured D —> K^K tv'^tv differential rate into components of definite T and CP transformation 
properties, projected onto the 0 angle and ttikk invariant mass. The uncertainties on the LHCb data points in Fig. 3(e-f) and 
Fig. 2(c-d) of Ref. [28] have been assumed equal to \/N -\- 8 and uncorrelated. 


sin ^ < 0. This allows for the marginalized differential 
distributions of definite T and CP properties to be de¬ 
rived. The left panel of Fig. 2 for instance shows those 
differential distributions projected onto the (j) angle (f.e., 
marginalized over the four other phase-space variables). 
The respective c^(j) and s2(j) dependences of the T-even- 
CP-even and T-odd-CP-even differential rates expected 
from a dominant contribution are clearly visible 
while the T-even-CP-odd and T-odd-CP-odd distribu¬ 
tions are roughly compatible with zero. 

An oscillatory pattern can however be distinguished 
in the T-odd-CP-odd differential rate. The A„ = 

/dd) signjsin n(/)} | asymmetries (see Table III) 

notably point at the presence of a sizable sin S4> con¬ 
tribution: the Ag departure from zero is of about 2.6 


standard deviations ( 2.0 standard deviations for A 2 and 
A 13 ). If genuine, this rapid oscillatory behavior would 
indicate the presence of a CP-violating phase difference 
but would not have contributed to the asymmetries that 
could be expected from a simple 4>p^ parametrization. 
Whether any resonance model considered as providing 
a fair description of that process would have included a 
contribution oscillating so rapidly is also unclear. 

E. Even more angular asymmetries 

Clearly, one way in which the presence of CP-violating 
phases could be probed without relying on a full descrip¬ 
tion of the dynamics of the process studied would be to 
evaluate systematically a wider range of triple-product 
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^ J 

1 27r 
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2\/2 

fdml 

dm^ 
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1 2-k 
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fdml 

dml 

J 

1 27r 

2ti 

2y/2 

fdml 

dml 

^ J 

1 27r 

2% 


fd$ signjcSa c9b c0} ^ 


fd$ sign{c6'a c9b sc;^} 


fd<l? sign{ c9a c9b ^ 


/ dr 

I d<l? sign{ c9a c9b 


f d<3? sign{ s20 } 


T-even 
CP-even 
T-odd 

d'J? CP-even 
T-even 
CP-odd 
t-odd 
CP-odd 
t-odd 


= + 


= + 


= + 


Af Oo a-j'i cos(So — Jjl) cos(v5o — (^l) 

i.j 

M a\ aj sin(5*L — ^o) cos((/3Y — tpi) 

i,j 

A/" ^ flo “'ll sin(5o - (5jj) sin(v5o - ^l) 

i.j 

M a\ aj cos(51 — 5^) sin(vpY — ^o) 
i.j 


d$ 


I d$ sign{ s2cf)} ^ 


I d<l> sign{ s20 } 


dr 

d$ 


fd<l> sign{c20} ^ 


CP-even 
t-even 
CP-odd 
t-odd 
CP-odd 
t-even 
CP-odd 

with Af = 


= + 


'dml 

dml 

2 -k 

2 k 

'dml 

dml 

2 k 

2 k 

'dm^ 

dml 

27 r 

2 k 

'dm^ 

dml 

27 r 

2 k 


Af aj_ a| cos(5V “ ) sin((^Y “ 

i,3 

^ X!“ii “ii 


2mo 87rmQ 


S-KrUa 


87rm^ 


TABLE II. Piecewise integrals from which information about the CP-conserving and CP-violating phases between different 
polarization amplitudes could be extracted, for a 0 —>■ (1 2) (3 4) decay involving spinless particles and proceeding through two 
intermediate vector resonances. 


n 


n 


n 


1 

-b58 ± 132 

6 

-bl64± 132 

11 

-fl28± 132 

2 

-259 ± 132 

7 

-flOl ± 132 

12 

-bl64 ± 132 

3 

-2 ± 132 

8 

-b337 ± 132 

13 

-b268 ± 132 

4 

-134 ± 132 

9 

-40 ± 132 

14 

-107 ± 132 

5 

-225 ± 132 

10 

-b41 ± 132 




TABLE III. A„ = fd$ sign{sinn(^} ^ P asymmetries 

t/ al* I CP-odd 

in the data collected by the LHCb Collaboration on the —>■ 

K~'K^'K~ decay. The uncertainties on the data points of 
Fig. 3(e-f) in Ref. [28] have been assumed equal to '/N + 8 
and uncorrelated. 


/o(ce) 

= 1, 

/i(ce) 

1 = c9, 

f2{c9) 

= 3c^6» - 1, 

Mc9) 

= c6i(3c^6» - 1), 

Mc9) 

= c6i(5c^6» -3), 

Mc9) 

= 5c^6' - 1, 

Mc9) 

> = 5c^9 -3, 

f7ic9) 

= c9 (5c^9 — 1), 

fs{c9) 

= c6i(3c^6» -1) (5c^9 -3), 

f9{c9) 

= (3c^9 - 1) (5c^9 - 1), 

fioic9) 

= c9(5c^9 -3) (5c^9 -1), 


asymmetries of the form 


TABLE IV. Natural set of functions of the 9a,b angles for 
the systematic construction of asymmetries in 0 —> (1 2) (3 4) 
decays involving spinless particles. 


I d$ sign{ fi{c9a) fm{ cOt) sin ncj)} 


dr 

d<i) 


t-odd 

CP-odd 


for all combinations of reasonably large integers /, m, and 
n. In the case of spinless final states forming two pairs of 
resonant intermediate states, the natural set of functions 
/ are products of the various c9 dependences arising in 



ax [GeV-i] 

<5x 

<px 

O 

II 

1 

1 

0 

II : 

2 

0 

0 

T : 

1 

1 

0.05 


TABLE V. Parameters chosen in the toy simulation of the 
D —^ (f)fP —>■ (A''^A'“)(7r^7r“) process. 
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spherical harmonics: 

1 = 1: c9, 

2: 3c26»-1, 

3: c9{5c‘^0-S), 5c‘^9-l, 

4: 35c46I -30c26»+3, c9 {7 c'^9 - S), 7 c^9 - 1, 

The dependences upon s9 have been dropped as they 
have no influence on the sign of the associated Legen¬ 
dre polynomials P™. Here, the set of / functions could 
therefore be defined as in Table IV, keeping in mind that 
another choice would be needed for final states carrying 
spin. 


F. Invariant mass dependence 



sin np sin rup 

sin 710 

(^TTTT ^p') 

sinm^ 

cos 9kk cos 9^^ 

n = 1 

-2.0 

0.86 

1.2 

-2.8 

2 

-4.0 

0.88 

3.4 

0.40 

3 

0.20 

0.15 

-0.014 

-1.4 

4 

0.30 

-0.014 

-0.52 

-1.5 

5 

0.30 

0.95 

-1.2 

-0.65 

6 

-0.40 

0.20 

1.0 

0.057 

7 

-2.0 

2.4 

2.4 

-0.042 

8 

-0.70 

0.37 

1.5 

0.27 

9 

-0.60 

-0.69 

0.88 

-0.75 

10 

0.60 

-2.2 

-0.78 

1.0 

11 

-2.0 

0.53 

2.1 

-0.15 

12 

-0.20 

-0.092 

0.55 

-1.7 

13 

-1.0 

0.30 

1.2 

0.67 

14 

0.20 

1.7 

0.30 

-0.70 


Let us still focus on the parametrization of the phase 
space privileging the 0 —> (12) (3 4) type of topology. 
Upon phase-space integration, the ^ invariant mass 
dependence of the decay rate could also lead to losses of 
sensitivity to CP-violating phases. This happens when it 
causes the T-odd-CP-odd piece of the differential decay 
rate to change sign. Guessing where this could happen is 
in general difflcult. However, when resonances are clearly 
identified, one at least knows the real parts of the asso¬ 
ciated propagators 


Re 


1 


ml-M^ + iTM } 


(m2-M2)2-Hr2M2 


change sign at the resonances (and could possibly appear 
in interferences). 

Once again, a glimpse at the LHCb data shows such 
a behavior actually occurs in the ttikk invariant mass 
spectrum, although in the T-odd-CP-even piece of the 
differential rate which is not directly relevant for the ex¬ 
traction of CP-violating phases (see right panel of Fig. 2). 

Therefore, when constructing asymmetries systemati¬ 
cally one may also wish to consider sign{m2 — Mf} and 
sign{m2 — Mj} as weight functions, for the known res¬ 
onances appearing at in the m^ ^ invariant mass 

spectra. 


G. Binned analyses 

Instead of constructing asymmetries, one may rather 
adopt the approach of Ref. [28] and bin the phase space. 
Care must however be taken in the binning choice. 
Putting together in one bin, regions of the phase space 
in which the T-odd-CP-odd part of the differential rate 
changes sign would result in sensitivity losses. 

These can be assessed using a toy simulation. We con¬ 
sidered massless kaons and pions and generated events 
using MadCraphS [56] with the following matrix elements 


TABLE VI. Departure from zero expressed in standard devi¬ 
ations for a few asymmetries, computed with the simulated 
sample of D —> (/)p° =>■ {K'^ decays. Only statis¬ 
tical uncertainties are accounted for. 


for the Dp(j), (j)KK, and pnir interactions: 


Dp(t) ■■ Cp (^Ao 


1 - 


2 2 
PpP^ 


{Pp‘P<p)‘^^ 




(j)KK : {pK+f, - Pk-^,) , 

pTTTT : ■ 


The linear polarization amplitudes described earlier are 
then 


Aq — 


_ Aq A(mj, _ 

12(m|.^ — m^+ im^V — nip -|- impVp) ’ 


All = 


—; 

KK 


-X(m\ 


- ^1 + im^T- mj, + impTp) ’ 


^ 6(m|f^ — mj, + im^jlT 4 ,) (m^„ — mj + impTp)' 

where each of the Aq y _l were given both a CP-even and 
CP-odd phase: Ax = ax for X = 0, ||,T. 

These parameters were fixed as in Table V and 40 000 
and decays generated. The decomposition of the 
(p differential distribution obtained is displayed in the left 
panel of Fig. 3. 

A larger magnitude for ay than for a± causes the T- 
even-CP-even piece of the differential rate to have a dip 
at 7r/2. The nonvanishing difference in CP-conserving 
phases (5y — (5j_ sources the sin 2(t) dependence of the T- 
odd-CP-even contribution. No structure is generated in 
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FIG. 3. Simulated D —>■ —>■ {K'^K~){'K'^n~) decay and partial rate decomposition in components of definite T and CP 

transformation properties as in Eq. (2), projected onto the (p angle and uikk invariant mass. Only statistical uncertainties are 
displayed. 


the T-even-CP-odd distribution while a small difference 
in CP-violating phases allows for a sin 2^ de¬ 

pendence in the T-odd-CP-odd differential rate. Due to 
limited statistics, the latter is barely visible in Fig. 3. Ad¬ 
ditionally, a sin 29a sin 20^ sin (j) dependence of each piece 
of the differential rate is washed out upon integration 
over the 9a,b angles. One can also notice a sign change 
in the T-odd-CP-even differential rate projected on the 
rriKK variable at the = 1.02 GeV resonance (see right 
panel of Fig. 3). 

Computing 


d$ sign fi ( c9a ) fm ( c9b ) sin ncf) 


W{ml - Mf) W{ml - M])\ 

7 -J 


dd) 


f-odd 

CP-odd 


( 6 ) 


asymmetries as prescribed earlier, one observes the ex¬ 
pected excesses for (Z, m; n) = (0,0; 2) and (1,1; 1). They 
are of 4.0 and 2.8 standard deviations, respectively (see 
Table VI, only statistical uncertainties have been ac¬ 
counted for). Using additional sign{m|f^ — to|} and 
sign{mj.„. — TOp} weight functions does not enhance the 
excesses’ significance. 

The LHCb Collaboration partitioned the phase space 
in 32 bins (two bins per kinematic variable) and esti¬ 
mated the combined departure from zero using a chi- 
squared test [28]. The separation between the two bins 
of the (j) variable was set at 1.99 rad (its domain is re¬ 
stricted to the [0,7r] interval here) and between —0.28 
and -1-0.28 for cos9kk and cos0,r7r- 

In our simulated sample, a chi-squared test with only 
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FIG. 4. The four components of the D —> K differential rate having definite T, E*, and CP transformation properties 

that could have been measured with an untagged sample. The uncertainties on the LHCb data points in from Fig. 3(a-d) of 
Ref. [28] have been assumed equal to '/N + 8 and uncorrelated. 


two bins of boundary located at 7r/2 in the (j) variable 
gives a departure from zero for the T-odd-CP-odd dif¬ 
ferential rate equivalent to 3.9 Gaussian standard devi¬ 
ations. Once the bins’ boundary is moved to 1.99 rad, 
this significance slightly diminishes to 3.7 tr. The loss of 
sensitivity is more significant for the only other binning 
relevant to this simplified simulation. With two bins in 
both the cos 9kk and cos directions, the T-odd-CP- 
odd differential rate departs from zero at the 2.8 cr level 
when the bin boundaries are chosen at 0, and at the 
1.1 CT level only when they are respectively taken at the 
extreme values of —0.28 and -1-0.28. 

The multiplication of unnecessary bins also leads to 
losses of sensitivity, in this scheme. With eight bins 
having boundaries at 7r/2 in the (j) angle and 0 in the 
cos 0KK,Tnr Variables, one for instance obtains an overall 


departure from zero of 3.5 standard deviations. 

H. Untagged D and B —> /f+JT“7r+7r“ samples 

Although a tagging of the has been carried out 
by the LHCb Collaboration in this —>■ 

decay, the self-conjugate final state could have motivated 
an untagged analysis. This is what was actually done in 
the study of the decay to the very same final state [8]. 

In both cases, the E* permutation defined in Sec¬ 
tion I.C sends {Ar+, AT”, 7r+, 7r“} to {AT”, Ar+, 7r“, 7r+}. 
In the parametrization of the phase space adopted thus 
far, it therefore acts trivially on the uikk, w,rir invari¬ 
ant masses, and on the (p azimuthal angle. The cosines 
of the polar angles in the K^K~ and tt+tt” subsystems 
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undergo the following transformations: 

E* [cos 9k+] = cos Ox- = — cos 6x+, 

E* [cos ] = cos 6^- = — cos . 

Practically, the untagged E*-odd distributions can there¬ 
fore be obtained by multiplying the weights of each 
recorded event by | sign{cos0;f+cos0,r+} and by con¬ 
sidering the absolute values of both cosines as kinematic 
variables. The same procedure carried out with the vari¬ 
able 4> yields the T-odd distributions. Using the LHCb 
measurement [28] , we display in Fig. 4 the projection onto 
the cos 9x+ and cos 0,r+ variables of the four differential 
rates that could have been measured with an untagged 
sample of D —>■ K'^ events. 

In its analysis of the 5° —>■ K'^ K~'K~^Tr~ decay [8], the 
LHCb Collaboration used a parametrization of the phase 
space privileging resonances in the and K~Tr'^ in¬ 

variant masses. A dominant K*^K* intermediate state 
motivated this choice. In that parametrization, the per¬ 
mutation E* exchanges the cosines of the polar angles 
defined in the two subsystems c0q and cOt , as well as 
their respective invariant masses nia and mb- Various E*- 
odd asymmetries can therefore be constructed by using 
weight functions g'(4)) (see Section I.D) proportional to 
either c9a — c9b , or rua — rrib- The E*-odd asymmetries 
measured in Ref. [8] were the ones possibly appearing for 
Kir subsystems forming partial waves of ja,b = 0 and 1. 
The arguments presented here to motivate the system¬ 
atic use of a wider range of T-odd-CP-odd asymmetries 
however also apply to T-even-E*-odd-CP-odd ones. 

III. CONCLUSIONS 

CP violation in K and B decays has so far been 
observed mostly through time-independent and time- 
integrated rate asymmetries. As multibody decays are 


being measured with an ever increasing accuracy, it is 
desirable to devote more attention to their rich differen¬ 
tial distributions. 

Taking, as an illustrative example, the —>■ 

K'^ decay whose differential distribution has 

recently been studied by the LHCb Collaboration [28], 
we proposed to measure a large set of generalized triple¬ 
product asymmetries. Their choice is guided by the 
topology—or resonance structure—of the contribution 
under scrutiny, by the spin of the particles involved, and 
by the location of the known resonances. An illustration 
of the procedure and of the losses of sensitivity that may 
occur with a suboptimal partition of the phase space was 
provided using a toy simulation. Such a procedure could 
obviously be applied to a wide range of other processes 
in which CP violation is searched for in differential dis¬ 
tributions. 

In charm decays, a signal of CP violation would clearly 
point at new physics. In B decays, however, standard- 
model CP violation is expected to be visible in some 
cases. We did not investigate whether cleaner probes for 
physics beyond the standard model could be constructed 
from differential observables. Clearly, more theoretical 
studies in this direction would be necessary. 

Our final point is to emphasize that more experimen¬ 
tal studies are needed in order to devise observables op¬ 
timized for specific processes. With the new data coming 
from LHCb and Belle H, such a task is timely. 
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